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Abstract 

By using type la supernovae (SNIa) to provide the luminosity distance (LD) directly, which 
depends on the value of the Hubble constant Hq = 100 h km s _1 Mpc -1 , and the angular diameter 
distance from galaxy clusters or baryon acoustic oscillations (BAOs) to give the derived LD ac¬ 
cording to the distance duality relation, we propose a model-independent method to determine h 
from the fact that different observations should give the same LD at a given redshift. Combining 
the Sloan Digital Sky Survey II (SDSS-II) SNIa from the MLCS2k2 light curve fit and galaxy 
cluster data, we find that at the lcr confidence level (CL), h = 0.5867 ± 0.0303 for the sample 
of the elliptical /? model for galaxy clusters, and h = 0.6199 ± 0.0293 for that of the spherical 
/3 model. The former is smaller than the values from other observations, whereas the latter is 
consistent with the Planck result at the 2a CL and agrees very well with the value reconstructed 
directly from the H(z) data. With the SDSS-II SNIa and BAO measurements, a tighter constraint, 
h = 0.6683 ± 0.0221, is obtained. For comparison, we also consider the Union 2.1 SNIa from the 
SALT2 light curve fitting. The results from the Union 2.1 SNIa are slightly larger than those from 
the SDSS-II SNIa, and the Union 2.1 SNIa + BAOs give the tightest value. We find that the values 
from SNIa + BAOs are quite consistent with those from the Planck and the BAOs, as well as the 
local measurement from Cepheids and very-low-redshift SNIa. 
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I. INTRODUCTION 


The determination of the Hubble constant Ho from astronomical observations is very 
important in cosmology, as it gives the present cosmic expansion speed and is related to 
the cosmic components and the size and age of our Universe. Its first accurate value, 
H 0 = 72 ± 8 km s -1 Mpc -1 at the lcr confidence level (CL) [1], was obtained using local 
measurements from the Hubble Space Telescope. About 10 years later, Riess et al. [2] 
improved this result significantly using Cephcids and very-low-redshift Type la supernovae 
(SNIa) and obtained H 0 = 73.8±2.4 km s' 1 Mpc -1 , where the uncertainty was reduced from 
~ 11% to ~ 3%. In addition to the cosmic distance ladder, the Hubble constant can also be 
determined using cosmic microwave background (CMB) measurement. On the basis of a six- 
parameter A cold dark matter (CDM) model, the nine-year Wilkinson Microwave Anisotropy 
Probe (WMAP9) data give H 0 = 70.0 ± 2.2 km s -1 Mpc -1 , a 3% determination [3, 4], 
Combining WMAP9 with smaller angular scale CMB data from the South Pole Telescope 
(SPT) and the Atacama Cosmology Telescope (ACT) as well as distance measurements 
derived from baryon acoustic oscillation (BAO) observations improves this determination to 
1.2% and gives H 0 = 68.76±0.84 km s -1 Mpc -1 [5], which is consistent with the result ( H 0 = 
68±2 km s -1 Mpc -1 ) determined by median statistics from Huchra’s final compilation [6, 7]. 
Recently, in the framework of a six-parameter ACDM model, the Planck Collaboration also 
found that H 0 = 67.3 A 1.2 km s -1 Mpc -1 [8] and H 0 = 67.8 ± 0.9 km s -1 Mpc -1 [9], which 
are consistent with the value (Ho = 67.3 ± 1.1 km s -1 Mpc -1 ) [10] obtained by combining 
galaxy BAO measurements from the Baryon Oscillation Spectroscopic Survey (BOSS) survey 
of Sloan Digital Sky Survey III (SDSS-III) [11] with high-precision measurements of the 
relative distances of the latest SNIa [12]. However, the Planck results are slightly smaller 
than the WMAP9 one and have a tension of about 2.5a with respect to local measurements 
given in [2], 

On one hand, this tension might indicate the need for new physics such as the cosmic 
variance [13] or a super-accelerating expansion [14], but on the other hand, it may also 
be due to systematic errors [15-19]. For example, revising the geometric maser distance 
to NGC 4258 from [20] and using this indicator to calibrate the data considered in [2], 
Efstathiou [15] yielded Ho = 70.6 ± 3.3 km s -1 Mpc -1 . Considering predominantly star¬ 
forming environments, Rigault et al. [16] found H 0 = 70.6 ± 2.6 km s -1 Mpc -1 . These 
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results are consistent with that of Planck. In addition, the systematic errors in the CMB 
analysis may also be responsible for some part of the tension. By removing the 217 x 
217 GHz detector set spectrum used in the Planck analysis, it was found that H 0 = 68.0 ± 
1.1 km s” 1 Mpc” 1 [19]. 

In addition to the very-low-redshift data from local measurements and the CMB data 
a,t z = 1089, the intermediate-redshift data provide a complementary tool for determin¬ 
ing the Hubble constant, and they help to reduce the uncertainties of H 0 significantly if 
they are combined with CMB data. Cheng and Huang [26] combined BAO data from the 
6dF Galaxy Survey [21], BOSS Data Release 11 (DR11) clustering of galaxies [11], Wig- 
gleZ [22], and z = 2.34 from the BOSS DR11 quasar Lyman-a forest lines [23], as well as 
simultaneous measurements of the two-dimensional two-point correlation function from the 
BOSS DR9 CMASS sample [24] and the two-dimensional matter power spectrum from the 
SDSS DR7 sample [25], and found that in a ACDM model, H 0 = 68.0 ± 1.1 km s” 1 Mpc” 1 , 
a 1.3% determination. Using the WMAP9+SPT+ACT+6dFGS+BOSS/DRll+i7o/Riess 
and basing their analysis on the six-parameter ACDM cosmology, Bennett et al. [5] obtained 
H 0 = 69.6 ± 0.7 km s” 1 Mpc” 1 , a 1% determination, which is the most accurate result to 
date. Through a nonparametric reconstruction of H(z) data, Busti et al. [27] obtained 
Ho = 64.9 ± 4.2 km s” 1 Mpc -1 by extrapolating the reconstruction to redshift 0. All these 
results seem to be consistent with that of Planck. However, using other intermediate-redshift 
data including angular diameter distances (ADDs) from the spherical f3 model galaxy cluster 
sample [28], 11 ages of old high-redshift galaxies [29, 30], 18 H(z) data points [31-33], and 
the BAO peak at z = 0.35 [34], Lima and Cunha [35] found Ho = 74.1 ± 2.2 km s” 1 Mpc” 1 
in a ACDM model. This value is consistent with the result obtained in [2], Later, Holanda 
et al. [36] considered the effects of different galaxy cluster samples and found that different 
samples give different results, i.e., H 0 = 70 ±4 km s” 1 Mpc” 1 if the spherical (3 model galaxy 
cluster sample [28] is replaced by the elliptical one [37]. Therefore, not only local and global 
measurements but also different intermediate-redshift data may yield controversial results 
for the value of the Hubble constant. 

The above discussion illustrates that, except for direct measurement from the cosmic 
distance ladder and nonparametric reconstruction from H{z) data, the reported results are 
based on the ACDM model. In this paper, we propose a new model-independent method to 
determine Ho by combining the observed luminosity distances (LDs) from SNIa data and 
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ADDs from galaxy clusters or BAOs. 


II. METHOD 

The measurement of SNIa is very important in modern cosmology, and it first indicates 
that the present expansion of our Universe is speeding up. The observed SNIa data are 
usually released in the form of the distance modulus /i, which is related to the LD D l by 

/i = 51 oSl 0 i^ + 25 - (1) 

The value of the distance modulus is determined mainly by the rest-frame peak magnitude 
m and the peak absolute magnitude M of SNIa. Because M is degenerate with Hq = 
lOOh km s -1 Mpc -1 , we must choose a value of h to obtain the distance modulus. For 
example, h = 0.65 in the SDSS-II sample [38], and h = 0.70 in the Union 2 one [39]. Because 
TA oc 7 in the Friedmann-Robertson Walker universe, letting h 0 denote the chosen value 
of h to determine the distance modulus, one finds that 

/i = - 5 log 10 , (2) 

U 0 

where /r 0 represents the released distance modulus with h = h 0 . Combining Eq. (2) with 
Eq. (1) gives the h-dependent LD from the distance modulus of SNIa: 

D L (h) = ^lO^ 5-5 Mpc = ^D lo . (3) 

Several light curve fitting methods are used to analyze SNIa data. Here, we consider the 
MLCS2k2 [40] fitting method and the SALT2 [41] method for comparison. In the MLCS2k2 
fitter, the estimated value and uncertainty of the distance modulus are obtained independent 
of any cosmological models. For this fitter, we use the 288 SDSS-II SNIa data points [38] in 
our analysis, which are given with h 0 = 0.65. 

For the SALT2 light curve fitting method, the released distance modulus is model- 
dependent, as the observed value of the distance modulus is obtained using fi = m — M + 
ax — /3c, where x is a stretch factor that describes the effects of the shapes of light curves on 
/i, and c is a color parameter that denotes the influences of the intrinsic color and reddening 
by dust, a, (3 , and M are fitted by minimizing the residuals in the Hubble diagram in the 
framework of the ACDM model. For the SALT2 fit, we consider the Union 2.1 [42] SNIa 
sample, which is obtained with h 0 = 0.7 and consists of 580 data points. 
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From Eq. (3), one can see that, if the LD can be obtained from other observations and is 
independent of h, the value of the Hubble constant can be determined directly by combining 
them with the SNIa observations, and the result is model-independent 1 . However, no 
observations other than SNIa can provide the observed value of the LD directly. Fortunately, 
the ADD Da can be obtained from some observations, such as galaxy clusters and BAOs. 
Using the distance duality relation between the LD and the ADD, D l = (1 + z) 2 Da, we 
can deduce the value of D l at a redshift of Da- Thus, if one has observed SNIa and galaxy 
cluster (or BAO) data at the same redshift and assumes that the distance duality relation 
is valid, the combination of these data provides a model-independent method to constrain 
h. This assumption is reasonable because the distance duality relation is valid in any metric 
theory of gravity, and its validity has been tested in [43-47] using different datasets. 

The above analysis shows that finding the observed LD (Dl 0 ) from the distance modulus 
of SNIa and the ADD from galaxy clusters or BAOs at the same redshift is necessary to 
constrain h in a model-independent way. Because the SNIa data sample is much larger 
than that of galaxy clusters or BAOs, to obtain a tight constraint, all ADD data should be 
kept. Owing to the usual absence of SNIa data at the same redshift as ADD data, we use 
a binning method to give the corresponding LD value of SNIa. This means that all SNIa 
data available in the range A z = jzp — za\ < A are binned, where za represents the redshift 
of the ADD data, and constant A represents the binned redshift region. Using the criterion 
that the standard deviation of the set of /io in a bin has the same order of magnitude as a ^ 0 , 
which is defined in Eq. (5), to determine the value of A, we choose A = 0.011. For all the 
selected data, an inverse variance weighted average is employed. If yUo, m denotes the rntli 
appropriate SNIa distance modulus data in a bin with the given redshift Zj, in which there 
are N t SNIa data points, and cr^ represents the corresponding observational uncertainty, 
one can straightforwardly find, using the conventional data reduction techniques in [48], that 



( 4 ) 




1 



( 5 ) 



where /xo,i represents the weighted mean distance modulus at a given redshift z 8 , and cr? 


1 Notice that for the Union 2.1 data, the result still depends on the ACDM model. 
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is its uncertainty. 

If the observational uncertainties of the distance modulus are released with the covariance 
matrix C, using a method similar to that given in [48], we find that the weighted mean 
distance modulus in a bin and its uncertainty can be expressed as follows: 

Era,n= 1 hO.mC 


hO ,i ~ 


J mn 


2 


Z—im,n= 1 mn 

1 


s-i 


-1 

mn 


( 6 ) 

(7) 

. In 


When only the statistical errors are considered, C is a diagonal matrix, and C nn = ct 2 q 
this case, Eqs. (6) and (7) reduce to Eqs. (4) and (5), respectively. 

Using Eq. (1), one can obtain the binned luminosity distance and the corresponding 
uncertainty: 


MQ,i 


D L0 ,i = 10 s’ 5 Mpc, 


( 8 ) 


CF D J ,n i ~ 


In 10 ^o.t . 

■10—-V,. 


•'LO.i 5 *” MO, 

Then, from Eq. (3) and the distance duality relation, the observed h 0 b s ,i can be obtained: 

h 0 D L0 , 


(9) 


hobs,i 


-(1 + Zi) 


-2 


Dx,i 

with °/L,i = h l D l% 1 + + h lDlo,i D Ki( l + y-yy,/ Here > is the observed 

ADD at the redshift and is its uncertainty. Minimizing the y 2 statistics 

( h ^o6s,z) 


( 10 ) 


X 2 = 


E 


erf 


( 11 ) 


we can obtain the constraint on h. The best-fit value is given by Xmim which is the minimum 
of x 2 . The la, 2 a, and 3 a CLs are determined through Ay 2 = y 2 — y< 1, Ay 2 < 4, 
and Ay 2 < 9, respectively. The probability density of h, P(h ) = Aexp(—y 2 /2), can be 
obtained, where A is a normalized coefficient. 

III. GALAXY CLUSTER SAMPLE 


In [49], it was suggested that the ADD of a galaxy cluster can be derived from a mea¬ 
surement of its intrinsic size on the basis of the Sunyaev-Zel’dovich (SZ) effect and X-ray 
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FIG. 1: The likelihood distributions of h from galaxy cluster + SDSS-II SNIa. The blue and purple 
solid lines represent the results from SDSS-II SNIa + 25 elliptical (3 model sample, and SDSS-II 
SNIa + 38 spherical f3 model sample, respectively. The shaded region shows the constraint from 
the latest Planck measurement with the la CL ( Hq = 67.8 ± 0.9 km s -1 Mpc -1 ) [9]. 

surface brightness observations. Here, we employ two different galaxy cluster samples: the 
spherical /3 model sample and elliptical /3 model sample. The former, compiled in [28], has 
38 ADD data points in the redshift region 0.14 < z < 0.89, which are obtained by using a 
hydrostatic equilibrium model to analyze the cluster plasma and dark matter distributions. 
The latter was compiled in [37] using an isothermal elliptical (3 model and consists of 25 
ADD data in the redshift region 0.02 < z < 0.79, in which 18 galaxy clusters were compiled 
in [49] and 7 galaxy clusters were compiled in [50]. 

A. Galaxy cluster sample ± SDSS-II SNIa 

The systematic errors in the SDSS-II SNIa data sample are not considered in our analysis. 
The binned distance modulus can be obtained from Eqs. (4) and (5). Fig. (1) shows the 
likelihood distributions of h, where the blue and purple solid lines represent the results from 
the SDSS-II SNIa + 25 elliptical (3 model sample and SDSS-II SNIa + 38 spherical (3 model 
sample, respectively, and the shaded region shows the constraint from the latest Planck 
measurement with the la CL ( H 0 = 67.8 ± 0.9 km s -1 Mpc -1 ) [9]. We find that at the la 
CL, h = 0.5867±0.0303 for the SDSS-II SNIa ± 25 galaxy clusters, and h = 0.6199 ±0.0293 
for the SDSS-II SNIa ± 38 galaxy clusters. 
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FIG. 2: The likelihood distributions of h from galaxy cluster + Union 2.1 SNIa. The blue and 
purple solid lines represent the results from Union 2.1 SNIa + 25 elliptical /3 model sample, and 
Union 2.1 SNIa + 38 spherical (3 model sample, respectively. For the SNIa data, only the statistical 
errors are considered. The shaded region shows the constraint from the latest Planck measurement 
with the la CL (Hq = 67.8 ± 0.9 km s” 1 Mpc” 1 ) [9]. 

B. Galaxy cluster sample + Union 2.1 SNIa 

In this case, we consider the Union 2.1 SNIa data without and with the systematic error. 
For the latter case, the covariance of the SNIa data in a bin is considered using Eqs. (6) 
and (7). Because the uncertainty of the galaxy cluster data is much larger than that of 
the SNIa data, the systematic error of the SNIa data has a negligible effect on the results. 
For example, at the la CL, h = 0.5940 ± 0.0303 and h = 0.5938 ± 0.0303 for 25 galaxy 
clusters + LInion 2.1 SNIa without and with systematic errors, and h = 0.6395 ± 0.0298 and 
h = 0.6393 ± 0.0298 for 38 galaxy clusters + Union 2.1 SNIa without and with systematic 
errors, respectively. Therefore, we plot only the likelihood distributions of h from galaxy 
clusters + Union 2.1 SNIa without the systematic errors, which are shown in Fig. (2), 
where the blue and purple solid lines represent the results from the Union 2.1 SNIa + 
elliptical (3 model sample and Union 2.1 SNIa + spherical /3 model sample, respectively. 
The shaded region shows the constraint from the latest Planck measurement at the la CL 
( H 0 = 67.8 ± 0.9 km s" 1 Mpc" 1 ) [9]. 

Comparing Figs. (1) and (2) reveals that, for the same galaxy cluster sample, different 
SNIa samples give almost the same result, but Union 2.1 favors a slightly larger h. The 
combinations of SNIa and the elliptical [3 model sample seem to favor a smaller Hubble 

8 









constant, which is marginally consistent with the latest Planck result at the 3a CL. The 
Hubble constants from SNIa + spherical (3 model are compatible with that from Planck 
at the 2 a CL and are quite consistent with that reconstructed from H(z) data directly 
given in [27], where H 0 = 64.9 ± 4.2 km s -1 Mpc -1 . It is easy to see that different galaxy 
cluster samples give different h constraints. This finding is the same as that obtained in 
[36] on the basis of a ACDM model in which galaxy cluster samples were combined with 
other data. However, our values for the Hubble constant are much less than those given in 
[36]. Thus, there is a slight incompatibility between the galaxy cluster + SNIa and Planck 
results, which is the same as that obtained in [8], where some inconsistencies between the 
cosmological parameter constraints derived from the Planck SZ cluster measurements and 
from the Planck CMB data were found. This suggests that assumptions about global cluster 
properties strongly affect the cosmological model inferred. 


IV. BARYON ACOUSTIC OSCILLATIONS 

BAOs (see [51] for a review) arise from the coupling of photons and baryons by Thomson 
scattering in the early universe. Competition between the radiation pressure and gravity 
gives rise to a system of standing sound waves within the plasma. At recombination, the in¬ 
teraction between photons and baryons ends abruptly because the free electrons are quickly 
captured, which leads to a slight overdensity of baryons at the sound horizon. The corre¬ 
sponding scale is the distance traveled by a sound wave in the plasma before recombination. 
This scale has been measured in the clustering distribution of galaxies today and can be 
used as a standard ruler. Combining measurements of the baryon acoustic peak and the 
Alcock-Paczynski distortion from galaxy clusters, the ADD data can be obtained. Note 
that, to obtain these data, a fiducial cosmological model is still needed. We summarize five 
low-redshift data points in Table (I) . Three of them are determined from the WiggleZ Dark 
Energy Survey [52], and the other two data points are from BOSS DR7 [53] and DR11 [54], 
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z Da(z) (Mpc) Survey 

0.44 1205 ± 114 WiggleZ[52] 

0.6 1380 ± 95 

0.73 1534 ± 107 

0.35 1050 ± 38 SDSS DR7[53] 

0.57 1380 ± 23 SDSS DR11 CMASS[54] 

TABLE I: Summary of the ADD measurements from the BAOs. 



FIG. 3: The likelihood distributions of h from BAOs + SDSS-II SNIa. The shaded region shows the 
constraint from the latest Planck measurement with the lcr CL (Hq = 67.8T0.9 km s _1 Mpc” 1 ) [9]. 


A. BAOs + SDSS-II SNIa 

The likelihood distribution of h from BAOs + SDSS-II SNIa is shown in Fig. (3). As in 
the previous section, only the statistical errors are considered for SDSS-II SNIa. At the la 
CL, we have h = 0.6683 ± 0.0221, which is tighter than the values from the galaxy cluster 
samples, although only five BAO data points are used. It is more consistent with the result 
from Planck than the values from the galaxy cluster samples. This consistency occurs at 
the lcr CL. 
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FIG. 4: The likelihood distributions of h from BAOs and Union 2.1 SNIa data. Blue and purple 
lines represent the results without and with the systematic errors of SNIa considered in a bin. 
The shaded region shows the constraint from the latest Planck measurement with the la CL 
(H 0 = 67.8 ± 0.9 km s~ x Mpc -1 ). 

B. BAOs + Union 2.1 SNIa 

Fig. (4) shows the likelihood distributions of h from BAOs and Union 2.1 SNIa. Blue and 
purple lines represent the results without and with the systematic errors of SNIa considered 
in a bin, respectively. A strong constraint on h is clearly obtained: h = 0.6905 ± 0.0176 
for BAOs + SDSS-II SNIa without the systematic errors, and h = 0.6803 ± 0.0187 for 
BAOs + SDSS-II SNIa with the systematic errors at the la CL, an approximately 2.5% 
determination. These results are consistent with those from Planck at the la CL, and the 
consistency improves when the systematic errors in a bin are included. 

The SDSS-II SNIa clearly favor a slightly smaller h, which is the same as that for the 
galaxy clusters. The values from SNIa + BAOs are larger than those from SNIa + galaxy 
cluster samples. They are quite consistent with the results from Planck with and without 
the systematic errors of the 217 x 217 GHz detector set spectrum considered [8, 9, 19] and 
the combination of all BAO data [26]. It is also compatible with the local measurement from 
Cepheids and very-low-redshift SNIa given in [15, 16]. 
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V. CONCLUSIONS 


In this paper, we propose a model-independent method to determine the value of the 
Hubble constant, which is a very important constant in cosmology. Because the distance 
modulus of SNIa depends on the chosen value of h, we obtain an h-dependent expression 
for the LD, which is given in Eq. (3). Thus, if the LD at a redshift can be determined 
by other observations, it can be combined with the one from SNIa to determine the value 
of h using the fact that two different observations should give the same value of the LD 
at a given redshift. According to the distance duality relation, the LD can be deduced 
from the observed ADD data. Using the SDSS-II SNIa released by the MLCS2k2 light 
curve fitter to provide the observed LD data and the galaxy cluster samples to provide the 
observed ADD data, we find that h = 0.5867 ± 0.0303 for the elliptical /3 model sample, 
and h = 0.6199 ± 0.0293 for the spherical /3 model sample. The former is smaller than the 
values from other observations, whereas the latter agrees very well with the value directly 
reconstructed from the H(z) data [27]. Thus, different galaxy cluster samples give different 
values of h. This result is the same as that obtained in [36] on the basis of a ACDM model 
with a combination of galaxy clusters and other observations. If the ADD data from BAO 
measurements are considered, we obtain h = 0.6683 ± 0.0221 at the la CL, which is larger 
than the values from SNIa + galaxy clusters. In addition, for a comparison, we consider 
the Union 2.1 SNIa sample given by the SALT2 light curve fitter. The results from Union 
2.1 SNIa are model-dependent, as the distance modulus from SALT2 is determined in the 
framework of the ACDM model. The Union 2.1 SNIa give a slightly larger value of h than 
the SDSS-II one. The covariance of the Union 2.1 data has an apparent effect on the result 
when BAO data are used. We find that the results from SNIa + BAOs are quite consistent 
with those from Planck [8, 9, 19] and the combination of all BAO data [26], as well as the 
local measurement of Cepheids and very-low-redshift SNIa [15, 16]. 
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